A Fast Transform for Spherical Harmonics
نویسنده
چکیده
Spherical Harmonics arise on the sphere S in the same way that the (Fourier) exponential functions {e}k∈Z arise on the circle. Spherical Harmonic series have many of the same wonderful properties as Fourier series, but have lacked one important thing: a numerically stable fast transform analogous to the Fast Fourier Transform. Without a fast transform, evaluating (or expanding in) Spherical Harmonic series on the computer is slow—for large computations prohibitively slow. This paper provides a fast transform. For a grid of O(N2) points on the sphere, a direct calculation has computational complexity O(N4), but a simple separation of variables and Fast Fourier Transform reduce it to O(N ) time. Here we present algorithms with times O(N 5/2 logN) and O(N2(logN)2). The problem quickly reduces to the fast application of matrices of Associated Legendre Functions of certain orders. The essential insight is that although these matrices are dense and oscillatory, locally they can be represented efficiently in trigonometric series.
منابع مشابه
A fast spherical harmonics transform algorithm
The spectral method with discrete spherical harmonics transform plays an important role in many applications. In spite of its advantages, the spherical harmonics transform has a drawback of high computational complexity, which is determined by that of the associated Legendre transform, and the direct computation requires time of O(N3) for cut-off frequency N . In this paper, we propose a fast a...
متن کاملFast spin ±2 spherical harmonics transforms and application in cosmology
An exact fast algorithm is developed for the spin-weighted spherical harmonics transforms of band-limited spin ±2 functions on the sphere. First, we recall the notion of spin functions on the sphere and their decomposition in an orthonormal basis of spin-weighted spherical harmonics. Second, we discuss the a priori O(L4) asymptotic complexity of the spin ±2 spherical harmonics transforms, where...
متن کاملFast Spin ±2 Spherical Harmonics Transforms
An exact fast algorithm is developed for the direct spin-weighted spherical harmonics transforms of bandlimited spin ±2 functions on the sphere. First, we define spin functions on the sphere and their decomposition in an orthonormal basis of spin-weighted spherical harmonics. Second, we discuss the a priori O(L4) asymptotic complexity of the spin ±2 spherical harmonics transforms, where 2L stan...
متن کاملA Fast Spherical Filter with Uniform Resolution
for transforming Legendre polynomial expansions, but it appears not to generalize to the spherical case. This paper introduces a fast algorithm for obtaining a uniform resolution representation of a function known at a latitude– Swarztrauber [12] has reviewed other transformation algolongitude grid on the surface of a sphere, equivalent to a triangular, rithms. isotropic truncation of the spher...
متن کاملUsing Fast Fourier Transform in the 3-d Multilevel Fast Multipole Algorithm
In this paper a method is presented how to perform interpolation and anterpolation in both spherical coordinates θ and φ by trigonometric polynomials and the fast Fourier transform (FFT) in the 3-D multilevel fast multipole algorithm (MLFMA). The proposed method is exact in interpolation and anterpolation, and has the high numerical efficiency of FFT. A numerical comparison suggests that the pr...
متن کامل